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1 Introduction 


Let y) be a measure (in general, complex), supported either on a finite set of products 
of curves in the complex x and y planes or, alternatively, on a domain in the complex 
plane, with the identifications {x = z,y = z). Consider the following 2A^-fold integral 

Ztv = y dy{xi,yi) ... j dy{xN,yN)‘^N{x)ANiy), (1.1) 

evaluated over some finite linear combination of support domains, where 

N N 

^n{x) = JJ(xi - Xj), A^iy) = JJ(|/i - yj) (1.2) 

i<j i<j 

are Vandermonde determinants. 

Such integrals arise, in particular, as partition functions and correlation functions for 
two-matrix models in those cases where one can reduce the integration over the matrix 
ensemble, e.g., via the Itzykson-Zuber, Harish-Chandra identity [14], to integrals over the 
eigenvalues {xi, of the two matrices. Depending on the specific choice of measure 

and support, these include models of normal matrices [5] with spectrum supported on 
some open region of the complex plane, coupled pairs of random hermitian matrices 
[14,20] or, more generally, normal matrices with spectrum supported on curve segments 
in the complex plane [3,4], and certain models of random unitary matrices [28,29]. (See 
Appendix^ for several examples.) 

From the viewpoint of perturbation theory, and also in order to apply methods from 
the theory of integrable systems to the study of such integrals, it is of interest to consider 
deformations of the measure of the general form 

d/i(x, y\t, n, m, t) := y), (1.3) 

where 

OO OO 

= a =1,2. (1.4) 

k=l k=l 

The four infinite sequences of complex numbers ...), ■■■)■, 

...), ...) and the integers n, m are viewed as independent 

deformation parameters. (Thus the quantities are not, in general, complex con¬ 
jugates of , t*^^)).) For brevity, we also use the notation t := t := 

The correspondingly deformed will be denoted Z 7 v(t,n,m, t). 
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In the present work we represent the integral (HID and its deformations as vacuum 
state expectation values (VEV) using two-component fermions, and use this to derive 
perturbation expansions for ZAr(t, n, m, f) as series in products of two or four Schur func¬ 
tions, with the four sets of continuous deformation parameters as their 

arguments. Such character expansions have been studied for various special cases of one 
and two matrix models by a number of authors [2,9,10,17,18,23,28,29] using a variety 
methods. The fermionic approach presented here is based on the original constructions 
of r-functions for integrable hierarchies (refs. [6,15] and [16]) as fermionic VEV’s . 

In section El after presenting the integral (HH) in the form of a fermionic vacuum 
expectation value, we show that the corresponding integral n, m, f) for the deformed 
measure (Q is a special case of a r-function for the coupled two-component KP hierarchy 
or, equivalently, a r-function for the two-component Toda lattice (TL) hierarchy [6,24-26]. 
(See ()2.15|1 below.) 

We also use this fermionic representation to derive the perturbation expansion for 
Zn (t,n, m,t) as weighted power series in the deformation parameters, written in the 
following alternative forms as sums over products of Schur functions 

Zjv(t,n,m,t) = ^ (1.5) 

= ( 1 . 6 ) 

A,/i 

= (1.7) 

A,/i 

= ( 1 . 8 ) 

A,/i 

= (1.9) 


where SA(t) is the Schur function corresponding to a partition A := (Ai,... A^v) of length 
N = i{X) (see [19], or Appendix iBl ) and the sums are over all quadruples (A,/i, i/, 77 ) 
or pairs (A, /i) of such partitions. The coefficients n, m), gxf^{N, n, m, t^^)), 

n, m, t(i),t(2)) and gxf!{N, are expressed as 

determinants of N x N submatrices of the matrix of bimoments 


Bik= / x"y'"dn{x,y), i,keZ 


and its deformations (see eqs. ()2.31|l - ()2.35jl and ()3.19j) below). 


( 1 . 10 ) 


3 




For normal matrix models, formula (USD was presented in [23]. The case of diagonal 
coefficients gx^ = Sx^rx was earlier considered in [9,10] in the context of models of pairs 
of hermitian matrices with Itzykson-Zuber coupling [14] and of normal matrices with 
axially symmetric interactions [5]. In the present work, we derive formulae (j1 .9j) 
by two different methods. First, in Section |21 via the two component fermionic calculus 
and then, in Section IHl by a direct calculation, as in [9,10,23], based on a standard 
determinant integral identity (the Andreief formula (13. Ill i. combined with the Cauchy- 
Littlewood identity (Q. We also present some of the earlier results in a more complete 
way; in particular, the expansions (HIEli-dllSI) are given for the case of coupled pairs of 
hermitian matrices. The fermionic representation to be used here however differs from 
the ones appearing previously in the context of matrix models in refs. [9,10,27], where the 
links with the one component KP (and TL) hierarchy were developed. It is also different 
from the fermionic approach of ref. [5], where two-dimensional fermions were used in the 
study of normal matrix models. 

Remark 1.1. One can obtain the A-fold integrals arising in one-matrix models as specializa¬ 
tions of the above, either by choosing the undeformed measure to contain a factor proportional 
to the Dirac delta function (5(x — y), in which case the matrix of bimoments becomes a Hankel 
matrix, or to (5(x — i), in which case it becomes a Toeplitz matrix. We shall not consider these 
specializations here. 

1.1 Free fermions 

Let A be the complex Clifford algebra over C generated by charged free fermions {/*, /ijigz, 
satisfying the anticommutation relations 


[U fj]+ = [fi, fj]+ = 0, [/„ /,]+ = (1.11) 

Any element of the linear part 

W := (©mezC/m) © (©mezC/m) (1-12) 

will be referred to as a free fermion. We also introduce the fermionic free helds 

fix) ■= fiv) ■= ^fky~''~\ (1-13) 

fcez fcez 

which may be viewed as generating functions for the fj, f/s. 
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This Clifford algebra has a standard Fock space representation defined as follows. 
Define the complementary, totally null (with respect to the underlying quadratic form) 


and mutually dual subspaces 

bFan := (©m<oC/m) © (© (©m>oC/m) © (©m<oC/m) , (1-14) 

and consider the left and right ^-modules 

F:=A/AWan. F-.= W,rA\A. (1.15) 

These are cyclic ^-modules generated by the vectors 

|0) = 1 mod AWan, (0| = 1 mod WcrA, (1-16) 

respectively, with the properties 

/m|0) = 0 (m < 0), fm\0) = 0 (m > 0), 

(0|/^ = 0 (m>0), (0|/^ = 0 (m<0). (1.17) 

The Fock spaces F and F are mutually dual, with the hermitian pairing defined via the 
linear form (0||0) on ^ called the vacuum expectation value. This is determined by 

(0|1|0) = 1; (0|/™/™|0) = 1, m<0; (0|/^/,„|0) = 1, m > 0, (1.18) 


(0|/n|0) = (0|/„|0) = (0|/™M0) = {0\fmfn\^) = 0; (0|/^/n|0) = 0, m^n, 

(1.19) 

together with the Wick theorem which implies, for any finite set of elements {wk G hF}, 

(0|wi ■ ■ -^271+110) = 0, 

(0|wi ■ • •w;2n|0) = ^ sgna{Q\wa{r)W^{2)\t)) ■ ■ ■ (t)\w„{2n-i)Wa{2n)\^)- (1-20) 

0'GS2n 

Here a runs over permutations for which cr(l) < cr(2),..., a{2n — 1) < cr(2n) and (j(l) < 
(j(3) < • ■ ■ < a{2n — 1). 

Now let {tai}i=i,...,Ar, be linear combinations of the //s only, j G Z, and 
linear combinations of the /j’s, j G Z. Then fjl.2()j) implies 

(0|wi ■ ■ ■WATtCAf ■ ■-wilO) = det ((OlwitCjlO)) |ij=i,...,Ar (1-21) 

Following refs. [6], [15], for all N E Z, we also introduce the states 

(iV| := (OICat (1.22) 
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where 



Cn 

_ f -pi^) 

— Jo ''' Jn -1 

if A > 0 

(1.23) 


Cn 

= f-i-■■ In 

if A < 0 

(1.24) 


Cn 

= 1 if A = 

0 

(1.25) 

and 







\N) :=Cm|0) 

(1.26) 

where 






Cn 

— In-I ■ ■ ■ fo 

if A > 0 

(1.27) 


Cn 

= In ■ ■ ■ I-1 

if A < 0 

(1.28) 


Cn 

= 1 if A = 

0 

(1.29) 


The states and (It:^ are referred to as the left and right charged vacnnm vectors, 

respectively, with charge N. 

In what follows we nse the notational convention 


AAr(a:) = det (xf (iV > 0), Ao(x) = 1, A 7 v(x) = 0 (iV < 0). (1.30) 

From the relations 

(0|/iv-fc/(x.)|0) = (0|/_^+fc_i/(7/,)|0) = 7/f-^ = 1, 2,..., (1.31) 

and (HHH), it follows that 

(7V|/(xi) • ■ ■ /(x„)|0) = <5„,;vAjv(a;), iV G Z, (1.32) 

(-7V|/(2/0---/(2/n)|0)=<5n,ivAAr(|/), AgZ. (1.33) 

1.2 Two-component fermions 

The 2-component fermion formalism is obtained by relabelling the above as follows. 

:= f2n+a-l , 

-\-oo 


■— fzn+a-l > 
-I-C50 

k=—oo 
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(1.34) 

(1.35) 







where a = 1,2. Then (mn) is equivalent to 

[/l”t /!'’']+ = l/i"t = 0. [/l"t /t®]+ = (1.36) 

We denote the right and left vacuum vectors respectively as 

|0,0):=|0), (0,0|:=(0|. (1.37) 


Relations (irmi then become, for a = 1, 2, 



o 

o' 

(m < 0), 

o 

IV 

o 

(1.38) 


o 

o 

(m > 0), 

o 

o 

A 

o 

(1.39) 

As in [6], [15], we also introduce the states 



( 


:= (0,0|C„(i)C'„(2), 

(1.40) 

where 


f (") 

) ■— Jo ■ 

r 

V 

o 

(1.41) 



f (") 

) •- J-i • 

A 

o 

(1.42) 



) := 1 

if = 0 

(1.43) 

and 


'nP‘^) 

:= C'„(2)C'„(i) 0, 0) 

(1.44) 

where 


._ f(«) 

^ if > 0 

(1.45) 



._ f(«) 

) • J n(°‘) 

• • /A if < 0 

(1.46) 



) := 1 

if = 0 

(1.47) 

The states (|1.4()|) and (11.4411 will be referred to. respectively as left and right charged 

vacuum 

vectors with charges 



It is 

easily verihed that 





/mV,*) = 0 

(m < n). 

V 

o 

IV 

(1.48) 


V*l/m^ = 0 

{m >n), 

V 

o 

iV 

* 

(1.49) 

and similarly 





= 0 

{m < n), 

* 

o 

IV 

(1.50) 


= 0 

{m >n), 

sly 

o 

A 

(1.51) 
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Remark 1.2. In subsequent calculations we use Wick’s theorem in the form (11.211) . In the 
two component setting we use formula (tr^ separately for each component. To calculate the 
vacuum expectation value of an operator O, first express it in the form 

o = Eo'‘'of. 


(1.52) 


Then 


(0,0|(9|0,0) = ^(0,0|(9f)of |0,0) = ^(0,0|0f^|0,0)(0,0|0f^|0,0) (1.53) 


(i)i 


l(2)| 


where Wick’s theorem in the form (tr^ may be applied to each factor (0, |0, 0). 

As a first application, note that Wick’s theorem and (usi-diisni) imply that 

k 

(A^,-A^l = (-l)^^(^+^)(5fc,7vA7v(a:)AAr(|/), (1.54) 

i=l 

where An{x) is dehned in fioni) . It then easily follows that 

k N 

{N + n, -lV-m| P[/<"(xi)/®( 9 j)|n, -m) = Ak);/) JJ x”(- 9 j)’”. 

(1.55) 


2=1 


2=1 


2 Fermionic representation for ^Ar(t, n, m, t) and double 
Schur function expansions 


2.1 Fermionic representation for Z^ 

Let 

41 :=/ f^^\x)P'>{y)dn{x,y). 

Then the following gives a fermionic VEV representation of Zn 

Zj^ = (-l)5^(^+i)(iv,-iV|A^|0,0) 
To see this, just use Oil) to write 


( 2 . 1 ) 

( 2 . 2 ) 


N 


(A,-A|A^|0,0) = (iV,-A| J] / f^^\x,)P\yMx,,y,)\0,0) (2.3) 

i=i 4 

r r ^ 

= (_l)|^W+i) ... ANix)ANiy)Yldy{xi,yi), (2.4) 


2 = 1 





where the last equality follows from (HMD. Now introduce 


OO A 1 . 

5 :^. (2.5) 

k=0 

Since 

{Ar,-Ar|yl^|0,0) = 0 if k^N (2.6) 

we may equivalently express (Q as 

= (-l)^^(^+i)iV!(iV,-iV|c/|0,0), iV>0 (2.7) 

By (HMD we also obtain 

{Ar,-Ar|e^|0,0) = 0, iV<0. (2.8) 

Next we show that for the deformed measure HI the expression (HID (or, equivalently 
(H2D) determines a r-function in the sense of integrable hierarchies. 

2.2 Z^(t,n, m,t) as a r-function [6] 

We begin by recalling briefly the important notion of r-functions as introduced by the 
Sato school. (See [6,15].) First, we dehne two inhnite linear families of operators H{t) 
H{t) by 

CXD OO 0000 


where the “commuting Hamiltonians” are bilinear combinations of fermion compo¬ 
nents of the form 


-hOO 


MO, 0 = 1 , 2 . 


( 2 . 10 ) 


Now, let 

g:=exp j ; : dp^^{x,y), 

a, 13 = 1 ,2 

where {dpay} is some 2x2 matrix of measures, and 


( 2 . 11 ) 


rp:=rp-<o,o\rp\o,o> 


( 2 . 12 ) 
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Then the expectation value 


rAr(t,n, m,t) = (A^ + n, —TV — —m), N,n,m&Z, (2-13) 

if it exists, is called the r-function of the two-component Toda Lattice (TL) hierarchy or, 
equivalently, the coupled two-component KP hierarchy. (See [6,26].) 

The sets of parameters t, t and also the integers N, n and m are called “higher times” 
of the two-component TL hierarchy. The parameters t = are higher times of 

the two-component KP hierarchy, and the second set t = is also a set of higher 

times, for a different two-component KP hierarchy. Together, they may be referred to 
as the coupled (two-component) KP hierarchy. The r-function (j2.13jl solves an inhnite 
number of bilinear (Hirota) equations. (Further details may be found in refs. [6,15,16]. 
We only note here that each equation contains a certain number of derivatives of the 
r-function with respect to the variables and a = 1,2, k = 1,2 ,..., and the 
equations may also relate r-functions with different values of the variables N,n,m.) 

To relate Z 7 v(t,n, m, t) to the r-function in the above sense, we choose the measures 
d/iii, dfi 22 and d^ 2 i to all vanish, and d^i 2 := dp. Thus g = of the form (ESD where 
A is defined by (12.11) . (In this case : A : just coincides with A.) We now will prove that, 
for At > 0, the resulting r-function, up to a simple explicit multiplicative factor, is equal 
to Z 7 v(t, n, m, t), where the measure dp is the deformed one in ()1.3|) : 

rAr(t, n, m, t) := {N + n, —N — —m) (2-14) 

= (2.15) 

N 

AAr(a;)A7v(p) JJ dn{xk, Pfc|t, n, m, f), 

fc=i 

where 

c(t,t):=e T . 

Proof: Eqns. f)B-llfl and ()B-12fl i of Appendix |B1 imply that 

,-am ^ ^v(..-\im)+v(r\ia>)fm(x)P\y) 

Using dehnition one can check that 

m) = |n,m), (n,= (n,m|, n,m 


(2.16) 

(2.17) 

( 2 . 18 ) 

(2.19) 

( 2 . 20 ) 


= (-l)i«(«+i)+>»Wc(t,t) 
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and, from the Heisenberg algebra relations 


= k5^^p5u,-u fc, / = ±1, ±2 ,..., «, /3 = 1,2, (2.21) 

it follows that 

gH(t)grr(t) _ ( 2 . 22 ) 

Combining these relations, we have 

(iV+n, —N—m\e^^^'^ge^^^^\n^ —m) = \n+ n, — A^——m), 

(2.23) 

where 

A{t,n,m,t) := J f^^'>{x)P'>{y)dij{x,y\t,n,m,t). (2.24) 

Finally, to remove the nnmbers n and m from the vacunm states we nse (|1-55D, giving 

{N + n,-N - -m) = e-S-iS"=i^4“^*t’(_i)-^(iV, 0) 

(2.25) 

Now (ITT3D follows from (IT^ in the same way that dZID was shown eqnivalent to (|22D. 


2.3 Perturbation series in the variables t(“\t(“) 

In this section we derive the expansions 


ZAr(t, n, m, t) = A^! 

5 ^ n, m, t( 2 ))sA(t^^^)s^(t(^)) 

(2.26) 


A,;/ 

£(A)/(^)<iV 


= N\ 

Y1 t( 2 ))sA(t(^))s^(t( 2 )) 

(2.27) 


a,m 

£(A),^(^)<iV 


= N\ 

Y1 t( 2 ))sA(t(^))s^(t( 2 )) 

(2.28) 


A,;/ 

£(A),^(^)<iV 


= N\ 


(2.29) 


a,m 

£(A)/(^)<iV 


where the snms range over all pairs of partitions A = (Ai,..., Aat), y = (pi,. 
lengths < N. For snch partitions, we also dehne the labels 

• •, I^n) of 

hi = 

Aj — i + A^, h[ = yi — i + N. 

(2.30) 
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In terms of these quantities, define four N x N determinants 



where c := (—and 



(2.35) 


is the matrix of deformed bimoments. 

The series show that we actually have a “quadruple” system of one- 

component TL hierarchies, since each of the double Schur series (I2I2EI)-(I2I2ni) has the 
form of a Takasaki expansion [24,25] of a one-component TL r-function, where the sets of 
higher time variables are, respectively, and One 

can also view Zjv(t, n, m, t) as determining a quadruple system of KP r-functions, in 
which each set, and plays the role of (one-component) KP higher times. 

We begin by proving First consider the case n = m = 0 and fd) = t^) = 0, 

then 


Zjv(t, 0,0,0) = -iV|e^d)^iV|g^g^_ 


(2.36) 


Rewrite 



(2.37) 


in terms of the component operators as 



(2.38) 


i,k£Z 


where 



(2.39) 


are the bimoments. 

Because of the anticommutation relations 


[/]■’./fl+ = o. Vj.iez 


(2.40) 
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and 


(2.41) 


(2.42) 


/ili|0,0) = 0, /f|0,0) = 0, fc>0, 

nothing is changed by making the snbstitntion 

A ^ / 1 ++ := ^ 

i,k>0 

(“projection to the creative corner”) when evalnating the r-fnnction. To compnte 

(2.43) 

we nse the relation 

(AT.• • vkhaT.io. °> = 

( 2 . 44 ) 

where A = (Ai,...,AAr) and p = (/xi,..., pat) are the partitions related to the labels 
{hi, hjj by (|2.3()jl . which is proved in Appendix IHl nsing Wick’s theorem. 

Note that each term (EH, for any given pair of decreasing sets of non-negative 
integers 

hi > h2 ■ ■ ■ > /lAT > 0, h{ > h2 ■ ■ ■ > h'j^ > 0 (2-45) 

occnrs in ()2.43|1 (A^!)^ times, mnltiplied by a product of the moments, with the sign 
determined by the permutations of the ordering, 

«»,(!),*(/.;) • ■ ■ ( 2 . 46 ) 

where a, a G are the permutations of the indices. (The index sets consist of distinct 
elements, since all VEV’s with repeated values of hi or h' vanish, and re-ordering distinct 
sets to satisfy changes the terms by the sign of the permutation.) Substitution of 

fl2.44|l into ()2.43|1 thus gives 

=iV! (2.47) 


X,fi 

£(X),£^l)<N 


where 


hv (^) = det iBh.^h>.), i,j = (2.48) 

for partitions with < N. Finally, we restore the dependence on the remaining 

variables (n, m, and higher times (t^^^t^^^)) by including it in the measure. We thus 
obtain 

{N + n, —N — A^(t) |n, —m) = N\ E t)sA(t^^^)s^(t^2)), (2.49) 

X,fi 

i(X),£{fi)<N 
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(2.50) 


where 2l(t) = y4(0,0,0,t) (see (12.241) ') and where 


0 + + 

yxfi 


(A^, n, m, = det (5 


hi+n,h'^+m 






if £(A),£(/i) < N. 

Relations are proved similarly. First, we eqnate the variables n, m and 

all irrelevant higher times to zero. (For the cases (j2.27j) . (j2.28j) and (j2.29j) . these times are 
(td),t(^)), and (tb),t(^)), respectively.) Then, for the cases (j2.27|) . (j2.28|) and 

dm, we consider 


and evalnate 


A - 

A__ 

. .(1) n(2) p 

(2.51) 



i,k>0 


A - 


f(i) A2) p 

•“ JN+iJk-N^N+i,l-k+N, 

(2.52) 



i>,k>0 


A - 

+ 

1 

t 

sr^ f(i) A2) p 

•— 2_^ /Ar_j_i/_fc_i_Ar-DW-i-l,fc+V 

(2.53) 


i>,k>0 


(2.54) 

(2.55) 

(2.56) 


For the evalnation of (EMD, (IT331) and (E3!n) we nse the following formnlae (see 
Appendix |Bj) 


(N, -A'l ■ ■ ■ A‘y„-i/4’-«e 

= (-l)"sA(E“')i>„(tW), 


h2) 


( 1 ) 


„iy(2)(t(2))_i^(l)(t(l)) 


| 0 , 0 ) 


(2.57) 

(2.58) 

(2.59) 


where A = (Ai,..., Aat) and p = (pi,..., ^n) are again the partitions related to the labels 
{hi, h'j} by ()2.30j) . 

Finally, we restore the dependence on the remaining variables {{n,m) and the appro¬ 
priate higher times, which are respectively and (4*-^^ t^^^) by inclnding 

these in the dehnition of the measnre, to obtain (IT7fl) . (IT^ and (EH) . 
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2.4 Generalized 2-matrix models with polynomial potentials 


In the case of the standard exponentially coupled hermitian two-matrix model [14], and 
a pair of polynomial potentials of degrees p -|- 1, g -|- 1, or the generalized case of pairs of 
normal matrices (2NMM) with spectra supported on some specihed curve segments (see 
refs. [3,4]), we may consider the unperturbed measure as corresponding to the leading 
monomial potentials 


dpi{x,y) := 


(2.60) 


/ p+i q+i r r 

{*)dii{x,y) = EE K.ab / / (*)e p+i ^"^^dxdy, (2.61) 

a=l b=l '^'"6 

with the integration contours {y^, Ffe} chosen in such a way that all the bimoment integrals 
are hnite: 

r r '1^ ^ 

x^y^dy{x,y) = EE i^ab / / x^y^e p+i 9+1 ^^^dxdy < oo. (2.62) 

i=l i = l ''"to- 

Using the notations of refs. [3,4], the deformations of the measure may be chosen to be 
exponentials of lower degree polynomials, defined by parameters ui,... ,Up and vi,... ,Vq, 
as well as multiplicative monomials in x and y of degrees n, m G N 

(i/i(a:, ?/|u, V, n, m) := (— l)'"a:"'p™'e“^^=i -fy'e~ p+i ~^^'^dxdy. (2.63) 

Formula (I2.26j) is applicable in this case if we put = 0 and dehne and 

as 


t(') = G := (-«i,...,0, 0, 0,...), 

2 p 

0 , 0 , 0 ,...). 

z q 

The partition function may therefore be expressed as 


72NMM 

-‘N 


(u, n, m, v) = N\ E gxf,{N, n, m)sx{tu)sp{K), 


(2.64) 

(2.65) 

( 2 . 66 ) 




where 


n, m) = det {Bh^+rL,h'.+m) ki=i,...,Ar, hi = \i - i + N, = yi - i + N, (2.67) 
and Bjk are given by ()2.62jl . 
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3 Direct derivation of Schur function expansions 

We now give an alternative derivation of formulae dUD-dlSl) through a direct reduction 
of the multiple double integral to determinants in the bimoments, using the well-known 
identity [1] 

N 

/ ya)det(j)i{xj)det'ijjkiyi), = A^ldetG, 1 < < N. (3.1) 

'' a=l 

where {0i,cire an arbitrary set of pairs of functions whose products are inte- 
grable with respect to the two variable measure dfi{x, y) on some suitable domain and 


Gij := J dn{x,y)(j)i{x)^pj{y) l<i,j <N. (3.2) 

We also make use of the following form of the Cauchy-Littlewood identity [19] 

^J2sx{t)sx{t), (3.3) 

A 

where the sum is over all partitions A. 

For any hnite set of variables {xi,... ,xn), we denote by [x] the inhnite sequence of 
monomial sums 


N 


N 


X := 




a=l 


a=l 


1 ^ 
a=l 


N 


xi,---), [a; ^] := E 


X 


-1 


a=l 


1 ^ 

- 

’ 2 “ ’ 
a=l 


1 ^ 
a a=l 


X 


-3 


(3.4) 

The identity (13.3j) may be used to express the exponential factors 
qV(x yt)^ ^ ^v{y yt( 2 )) series in Schur functions: 




gV(x-yt(i)) 


A 

^SA([a:"^])sA(t(^)), 


A 


(3.5) 

(3.6) 


where the sums are over partitions A with length ^{\) < N. 

Substituting the expansions (j3.5|) into the matrix integral dehning Zn (t, n, m, t) gives 


N 


a=l 


Z^r(t, n, m, n / dy,{xa,ya\t,n,m,t) I AN{x)AN{y) 
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N 


a=l 


JJ ( / dlJ,{Xa,ya\0,n,m,t)\AN{x)AN{y)^Sx{t^^^)sx{[x])^Sf,{t^^^)Sf^ 


Making use of the Jacobi-Trudy formula 

\j-j+N. 


where 


gives 


where 


A„(i) ’ '‘'‘I' A„(i) ’ 

A := Ai > A 2 > ■ ■ ■ > Ai(A), /i := /ii > /i2 > ■ • ■ > 

A 11 


(3.7) 

(3.8) 

(3.9) 
(3.10) 


7”“(t) ■= II ( det(i 

a=l A-f / 


i^-''+'')det(!,«-+''). 


/i./—Z+A /’n 


(3.11) 


Applying the identity (EH) for := x\ (j)k{y) ■= y^ then gives 

47 = iV!detG-, (3.12) 

where 

(G'v )*i := j (0, 0, (3.13) 

proving the equality (Eini). The other three cases of eqs. are derived simi¬ 

larly. 

We now derive the quadruple Schur function expansion (EH). To do this, we use the 
product formula for characters (see [19]) 


Sa(W)s,(H) = 5^Cvi>o(W). 


(3.14) 


where the the Littlewood-Richardson coefficients {c"^} are the multiplicities with which 
the tensor representations of type a appear in the tensor product of those of types A and 
y. We have, as a consequence of the Jacobi-Trudy formula. 


^]) = (xi • •-a^jv) '"lS£;([x]), 
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s,(l!/ ‘D = (!/i'"9w) '“Sii(l!/|) 


( 3 . 15 ) 


where v = (i/i,..., v^q) and r] = (r/i,..., //tv) and 

'■= — ^N-i+lt Vi ■= V'l ~ VN-i+l- (3.16) 

Combining this with 

SA(N)sp([a;]) = 5^c^ps«([x]), s^([|/])s^j([y]) = '^c^^fiS0{[y]), (3.17) 

a f} 

gives 


Z 7 v(t, n, m, t) 


(3.18) 




where 

Ixt,un{N,n,m) := ^ 

(3.19) 


a,l3 


and 

:= iV!det(5)^+„,K+^)|i<,,,<^, 

(3.20) 


li := tti — i + N — Ui, 1'- := Pi — i + N — rji. 

(3.21) 


4 Summary and further related work 

We have given a new two component fermionic representation for 2iV-fold integrals of 
type (HU and nsed it both to show their relationship to the two-component TL hierar¬ 
chy, and to get donble and qnadrnple Schnr fnnction expansions as pertnrbation series 
in the deformation parameters, generalizing analogons formnlae earlier obtained in [9,10] 
and [23]. These resnlts were also derived throngh a “direct” method, based on standard 
determinantal identities and character formnlae. In another work [11] we show how to 
nse this fermionic representation to rednce the evalnation of mnltiple integrals of rational 
symmetric fnnctions arising as determinantal correlation fnnctions in matrix models, ns- 
ing Wick’s theorem, to determinantal expressions involving at most donble integrals, and 
the associated biorthogonal polynomials. These resnlts may also alternatively be derived 
throngh direct methods [12], based on standard determinantal identities and partial frac¬ 
tion expansions. In [13] a fermionic representation for models of matrices conpled in a 
chain [7] is also derived. 
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A Appendix: Relation to unitary, hermitian and normal 
matrix models 

In this appendix, we give a number of examples of matrix integrals that reduce to 2N fold 
integrals over their eigenvalues of the form (EH). Two-matrix models involve integrals of 
the type 

<F> = -^ [ F{Mi,M2)dn{Mi,M2), 

J 

Zn = j dn{Mi,M2) (A-1) 

where F is some function of the entries of the matrices Mi and M 2 and dfl{Mi, M 2 ) is some 
measure (possibly complex) on the set of such pairs of matrices. From among the various 
models that have been studied, we review here three types: (A) normal matrix models, 
with spectrum supported in open regions of the complex plane; (B) models of pairs of 
hermitian matrices or, more generally normal matrices with spectrum supported on curve 
segments in the complex plane and (C) models of unitary matrices. The problem is to 
reduce the ~ integrations over the independent matrix entries to just 2N integrations 
over the eigenvalues of Mi and M 2 . This is possible for certain choices of matrix measures. 

(A) The simplest are models of normal N x N matrices M, where 

M+ = Ml =M, [M, M+] = 0, (A-2) 

and the integrand is invariant under U(A^) conjugations. One can then diagonalize the 
matrix using such transformations, 

M^UMU\ U e V{N) (A-3) 

and integrate over the group U(A^) to reduce the integral to a multiple integral over the 
eigenvalues with a suitable domain of integration in the complex plane, and the induced 
product measure multiplied by the factor |AAr( 2 ;)p . The resulting reduced measure ()1.3|) 
can be a rather general one in the complex plane. 

For illustrative purposes, we consider the special case when the undeformed measure 
depends only on the combination MM^, and the deformed partition function is of the 
form 

Zm = j dM j dM'^(detM)”(detMt)'"e^^^^^^^^+^*-i^“^*+*»^'’^^^^‘). (A-4) 
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Diagonalizing by U(A^) conjugation and integrating over the group U(iV), this reduces, 
up to a proportionality constant Cat, to: 


N 


Z]\[ — Ctv / dzidzi 






(A-5) 


a=l 


Because of the polar rotational invariance in each complex ; 2 a-plane, the angular parts 
of the bimoment integrals appearing in (nmn) may be evaluated. Assuming, e.g. that 
n > m, the only nonvanishing terms in the sum (ITTUD are those for which the partitions 
A and /i are related by; 

Xi + n = fii + m (A-6) 

Evaluating the angular parts of the bimoment integrals, (IXTUl) reduces to: 

Zn = Cn 5fA(n)sA(t^^^)sA+n-m(t^^^) (A-7) 

A 


where A + n — m denotes the partition (Ai + n — m, ••• , Xj^ + n — m) and 


9\{n) 


.= TT 



gV(X)^y-i+n+7V^^ 


(A-8) 


(B) In the case where {Mi, M 2 ) are a pair of independent hermitian matrices, or normal 
matrices with spectral supports along some specihed curve segments, the problem is more 
involved. The hrst example is the hermitian two-matrix model of Itzykson and Zuber [14] 
(see also [20], where these are referred to as unitary ensembles). The partition function 
is 


In — 


gTr(yi {M,)+V2 (M 2 )) gTrMi M 2 dVt{Mi)dn{M2) 


(A-9) 


where 


N N N 

dQ{M) = JJ dMii JJ d’ltMij JJ d^Mij. 

1=1 i<j i<j 


(A-IO) 


The potentials Vi,V 2 can be fairly general, but most often are taken as polynomials. 
Diagonalizing the matrices Mi and M 2 via two distinct U(A) conjugations. 


Mj = DjXf/j \ i = 1,2, X := diag(xi,... ,a;Ar), F := diag(?/i,..., j/at), (A-11) 
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and integrating over Ui^ one obtains 


N 


In = V^ I I dx.dyi--- I / 

(A-12) 


i=l 


where 


Jn{X,Y) := 


e^rUXWY^jj 


'U(iV) 


(A-13) 

is the remaining integral over the nnitary gronp, dJJ is the normalized Haar measnre and 

(2^)iv(7v-i) 


Vn-.= 


nf=i k\ 


2 • 


(A-14) 


The observation of [14] was that 

f .TrVXWY^^fj = m 


^N-1 


det(e^* 


Vj' 


'U(Ar) 


. fc=l 


AAr(x)Ajv(?/)‘ 


(A-15) 


Using the anti-symmetry of determinants and changes of variables, one then obtains the 
reduced integral 


In — 


nf=i k\ 



N 


dx2dy,... / (A-16) 


i=l 


which is proportional to (HID for measures m of the form 

dfx{x,y) = dxdy. 


(A-17) 


This same computation is valid for any family of matrices that are unitarily diagonal- 
izable, even if the spectral support is not on the real axis. Thus, if (Mi, M 2 ) are taken 
as normal matrices with spectrum supported on some union of curves {7a}, for Mi and 
{Tb} for M 2 , then the reduced form of the partition function is obtained by replacing 
the double integrals over the real axes in ()A-12jl by integrals of the form given, e.g., in 
eq. (imi) for the case of polynomial potentials ( [3,4]). 

Generalizations of this construction were considered in [28] and in [22,23]. We may 
replace the interaction term dMD by a more general one of the form 

Tr{N,M,M,):= d^,Nrx{N)s^{M,M,), (A-18) 

A 

i{X)<N 
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where {r(j)}jgN+ is some sequence of complex numbers, 

^a(A^) := n r{N + j-t), (A-19) 

ij'eA 

with the product ranging over all nodes of the Young diagram A with positions and 


dx^N — sa(Iv) (A-20) 

is the dimension of the representation of GL(Y) given by irreducible tensors of type A. 
(We use here the abbreviated notation Tr{N, M 1 M 2 ) to express what, in the more general 
setting of [9,22,23] was denoted Tr{N,lp^, M 1 M 2 ). Note that expressions such as (IA-1811 . 
viewed as functions of the trace invariants of the matrix Mi M 2 , are also KP r-functions.) 
Here the Schur functions sx{Z) are interpreted as conjugation invariant functions dehned 
on the space of complex N x N matrices Z (see Appendix^l)- Then Jn{X, Y) is replaced 
in eq. (jA-12jl by the following more general integral (see [22,23] for details). 


Jiv.r-(X,Y) := / Tr{N,UXU-^Y)dM = V rASA(X)sA(H) 

JV{N) V 

^(A)<iV 

_ det(r^(l,Xi?/j)) 

“ A^{x)A^{y) ’ 


(A-21) 

(A-22) 


where 


00 


Tr{l, x) = 1 + ^ r(l) ■ ■ ■ r{k)x’^ 

(A-23) 

k=l 


c ■ ^ 

• T-rV-l T-rfc / -N • 

nfc=l 11^ = 1 ^0) 

(A-24) 


The reduced integral is again of the form (EH) with the measure ()A-17jl replaced by 

dy{x,y) = e^^^^'^^^'^^'^\r{^,xy)dxdy. (A-25) 


Note that, setting Y = Iat on the left hand side of eq. (EHH), and taking the limits 
{yj^l}j=i,...,N on the right gives 


rr{N,X) 


CN,r 

Uti k\ 


(A-26) 
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Choosing 


we 


obtain from (jA-26jl 


A,x)=e^, r{j) =- 
J 

(A-27) 

Tr{N,X) = ^dX\ 

(A-28) 


which is the Itzykson-Zuber case. 

Another example is obtained by choosing 

1 

Tr{l,x) = 


(1 — ’ 

for some pair of constants (a, z). In this case 

Uti k\ 


r{i) = 


z{a-N + j) 


C]\f^r 1 




where 


(a-7V + l),:= J](a-iV + j) 

j=i 

is the Pochhammer symbol, and 


JnAX,Y) = 


det(l — zXiUj) 


N-a-l 


/S.N{x)AN{y) 

Eq. ()A-26jl therefore gives 

rr{N,X) = det(lAr - zX)-A 
Thus the integral (IA-12|1 is replaced by 


(A-29) 


(A-30) 


(A-31) 


(A-32) 


(A-33) 


1^,. = j y e'^"(^i(^i)+^2(^2»det(l7v-zMiM2)-“dff(Mi)dfi(M2) 

' 27 r 2 \i^(^-b 1 r f JL pVi{xi)+V2{yi) 




J J A^(a:)A^(|/) 


zXiHi 


\a-N+l 


dx^dyii 

(A-34) 


which is proportional to (HH) for measures of the form 

dAx, y) = (1 - zxyf-'^-Axdy. (A-35) 
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Remark A.l. For the case a = 0 the right hand side of HA-821) . up to the constant factor 
is 

ANix)AN{y) 

which itself is a constant in the variables {xi,yi}, as of course it must be, since there is no 
coupling in this case between the pairs of matrices (Mi, M 2 ). This may be seen by noting that 
the numerator and denominator are polynomials of the same degree A — 1 in all the variables, 
and antisymmetric under the interchange of any pair Xj Xj or ?/j t/j. This is therefore a 
rational function without poles, hence a polynomial, which is bounded in all the variables, and 
therefore a constant. 

(C) For unitary two-matrix models [23] (referred to in [20] as circular ensembles). Mi G 
U(A), i = 1, 2 , we have 




jgen _ 


'V{N) J\J{N) 


grr(yi(Mi)+y2(M2))^^('^^ (MiM2)"^)d*Mid*M2 


i.dxi dyi 


CN,rVNN\ / ... /Aiv(a;)Aiv(l/)ne'"i("^)+'"^("“V,(l,(a;,2/,)-^)^^(A-37) 
J J Xi yi 


where d^Mi, d^M 2 are the Haar measure on two copies of the group U(A), and Xi and 
yi, i = 1,... ,N are eigenvalues of the matrices Mi and M 2 , with values on the unit circle 
in the complex Xi and yi planes. This is related to measures in m of the form 


dy = {xy) ^) 


dxdy 

xy 


An example is the case 


r,(A, (MiMa)-^) = 

a(1, {xy)-^) = 


(A-38) 


(A-39) 

(A-40) 


which was considered in detail in [28,29], where the large N limit of this model was also 
studied. 

Another simple example, which is the unitary analog of the one above, is 


jgen _ 



gTr(yi(Mi)+y2(M2))^g^j^ - z{MiM2)-^)-M,Mid,M2 
= CN,rVNN\ 


/ \ 1 ^ e^d^d+V2{yi) dxidyi 
Ajv(x)Ajv(y) n n - > (A-41) 


J- (1-^Ua-AT+l) ^ ^ ’ 

=1 V ^ 
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which corresponds to measures in m of the form 


dll = e^i(^)+^hy) 


xy J xy 


In this case, r is the same as (IX^ 


z{a-N + ]) 

r{j) := -^- 

J 


and hence 


—1\—a 


Tr{N, (M 1 M 2 )-') = det(I^ - z{MiM 2 )-^) 


Tr{l, —) = 

xy 


n _^')a-Af+l ■ 

' xy ^ 


(A-42) 


(A-43) 


(A-44) 

(A-45) 


B Appendix: derivation of formula (l2.44l) 


In this appendix, we recall some dehnitions about Schur function (see [19] for further 
details) and results from refs. [6,15] that are needed in the derivation. 

A partition is a sequence of non-negative integers in weakly decreasing order: 


A — (Ai, A 2 ,..., Ar,...) , Ai > A 2 > • • ■ > A,. > ..., (B-1) 

where we identify (Ai,..., A„) with (Ai,..., A„,, 0) and, in the case of an inhnite sequence, 
assume A^ = 0 for r 0. The Aj in (lED are called the parts of A. The number of 
nonzero parts is the length of A, denoted by ^(A). We use the notation t = {ti, ^ 2 , • • •) and 
t = {ti,i 2 ,...) in what follows. The elementary Schur functions {si(t)}jgN are dehned 
by 

exp(I/(x, t)) = ^ Si(t)x^, (B-2) 

k>0 

where 

00 

V{x,t) (B-3) 

k=l 

The Schur function SA(t) corresponding to a partition A is given by 

SA(t) = det(sA,-i+i(t)),<. , (B-4) 
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where, for A; < 0, we put = 0 . When the sequence is identihed with the 

elementary trace invariants {^Tr(Z'^)}jgN of elements Z G GL{n) or Z G U{n), the Schur 
function SA(t) is the trace (character) of the rank i{X) irreducible tensor representation 
whose symmetries are given by the Young diagram associated to the partition A. When 
interpreted in this way, as a class function on the group GL{n), its values are denoted 
s\{Z), and this may be extended, by continuity, to all complex N x N matrices Z. 

Using the notations of subsections 11.11 and 11.21 we dehne, for m ^ 0, 

(B-5) 

iez 


It follows that 

and also that 

and hence 


For any a G Al, let 


[Hmi fi\ fi—m) [Hmi fi] fi+my 


Hm\0) =0, m > 0, 

^ |g^_ 


(B-6) 

(B-7) 

(B-8) 


a(t) := ^ exp(ad^ Hmtm) a, (B-9) 

m=l 

oo 

a(t) := = exp (ad a- (B-10) 

m=l 

Using (jR-fij) it is easily verihed that 

/(a:)(t) = exp(U(a;,t))/(a;), f{y)it)=exp{-V{y,t))f{y), (B-11) 

/(x)(t) = exp(U(a;"\t))/(x), f{y){l)=exp{-V{y-\l))f{y), (B-12) 


and hence 


/i(t) = X]sfe(t)/i_fc, 

fc >0 

/i(t) = 5Zsfc(t)/i+fc, 

fc >0 

Dehne the states 


/»(*) = '^^k{-'t)fi+k, 

k>0 

/i(t) = 5Zsfc(-t)/i_fc. 
fc >0 


(B-13) 

(B-14) 


{iV| := {0|/„ ■ ■ ■ {-iV| := {0|/_, ■ ■ ■ 


(B-15) 
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Now consider the integers hi > ■ ■ ■ > > 0 related to a partition A = (Ai,..., Aat) by 


hi — Xi — i N. 


(B-16) 


It follows from ()B-8(1 that 

""‘"A. • ■ ■ AJO) = {Af|A.(t) ■ ■ ■ A„(t)|0). (B-17) 

By Wick’s theorem and eqs. ()B-13j) . (jB-14jl this is equal to 

det {{0\fN-ifhj{'t)\0)) \i,j=i,...,N = det SA^._j+i(t) |jj=i,...,Ar = SA(t), (B-18) 

where the second equality follows from (ira . Thus we obtain 

■ ■-/h^lO) =SA(t). (B-19) 


Similarly, we have 

■ ■ -/-A^-ilO) = det (0|/i_jv-i/-A,-i(t)|0) |*j=i,...,7v 
det SAj— t) 

= 5A(-t). (B-20) 

To prove eq. (ICTl we use factorization, (IbTI and (IB^ 

{N, _ /d)y42) ^ 

= ^ ■ ■ ■ f^\_,\0) 

(B-21) 

= (_l)|A^(^+i)s,(t(i))s^(t(2)), (B-22) 


where 

Similarly, from 
it follows that 

We therefore have 


hi := Xi — i + iV, h[ -.= Hi — i + N. 


{±N\H_^ = 0 

= (±A^| 


(B-23) 

(B-24) 

(B-25) 


(^l/Af-Ajv-l ■ ■ ■/a^-Ai-IC — (A^|/Ar_/ijv-l(t) ■ ■ ■/v-Ai-l(t)|0) (B-26) 
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By Wick’s theorem this is equal to 


det ((0|/i_i/jv_/,^._i(t)|0)) \i,j=i,...,N = det |ij=i,...,7v = SA(t), (B-27) 

where the last equality follows from (IB-14f) . Thus 

(B-28) 

Similarly, we obtain 

{ — N\fhN-N---fhi-Ne = det (0|/_i//i_Ar(t)|0) |jj=l,...,Ar 

det SAj— j+i( f) \i,j=l,...,N Sa( t), 

(B-29) 

and hence 

{-N\U,-N---fh^-Ne^-^^-^-\0) = (B-30) 

By factorization of the appropriate VEV of two-component fermions into the product 
of VEV’s of one component fermions, eqns. (IB-19|1 . (IB-20|1 . (IB-28|1 and ()B-30|1 then imply 
(EH, as well as (0371) . (173^ and (I73ni) . 

Remark B.l. Note that all partitions A considered in this appendix have length i(X) < N. 

C Appendix: Quadruple Schur function series for 
T-functions of two component Toda lattice and 
derivation of formula 

Using results from refs. [6,15] one proves (see e.g. [9]) the formulae 

(Ar|eS^=i^'»*'= = 5^SA(t)(iV;A|, 5A(t)|iV; A), (C-1) 

A A 

where the sums range over all partitions and the Fock vectors on the right hand sides are 
dehned as follows 


(7V;A| :=a(iV,A)(iV| 

i=l 

(C-2) 

k 

|iV;A) :=a(iV,A)n/iV-ft-i/^+aJ^^)- 

(C-3) 


i=l 
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where 


a{N, A) = (_i)/3i+-+/3^=+5^(^-i) (C-4) 

Here k is the number of diagonal nodes of the Young diagram corresponding to A, aj is 
the number of nodes in the jth row to the right of the {j,j) diagonal one and f5j is the 
number of nodes in the column below it. It follows that the sets of integers {ai}, {/3i} are 
strictly decreasing (ai > • • • > and (/3i > ■ ■ • > j3k), and they uniquely determine A. 
The partition A is expressed in Frobenius notation (see ref. [19]) as: 

\ = {ai,... ,ak\Pi,..., f3k). (C-5) 

The transposed partition is denoted A*’’, which in Frobenius notation is 


A**^ {Pi, .. ., Pk\oii, ..., ak). 


(C-6) 


From the Jacobi-Trudy formulae (UHl) it follows that 

f{xi)---f{xN)\0) = ^ \N-,X)sx{[x])An{x) (C-7) 

A 

i{X)<N 

f{yi)---f{yN)\o)= I-^;^)(-i)'^'sAf([i/])A7v(2/). (C-8) 

A 

e{X)<N 

In the two component setting, we use the following notations: 




fc(2) 


fc(i) 


:= a(Af<‘>, A<‘>)a(Af®, A®) \ J] 


fc(2) 

ilf/J 




2 = 1 


fc(i) 


2 = 1 


2 = 1 


NW+ai^NW-Pi-V 

(C-9) 


:= a(iV'y \<^>)a{NP->, A'«) J] 11 1 


fc(2) 

(1) ,(i) TTf(2) 


NW-Pi-l-’ NW+Oi 11 Ar(2)-/3i-i^Ar(2)+, 

2=1 2=1 


(C-10) 


where 

Ad) = {qA\ ..., \j3[^\ ..., j = 1,2. (C-11) 

The vectors {{N^^\ {| A^^\ form dual orthonormal 

bases for the left and right Fock spaces: 
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Using ED for each component, we then have 

+ A.iVPI +m|, 

x,fi 

(C-13) 

eEf..«iX’-E£,Ha.i‘>|iv<‘),]vM) = +A,]VPI +^). 

A,/i 

(C-14) 

This implies that each two-component TL r-function (imn) may be expanded in a quadruple 
Schur function series: 

{N+n, ge^^^'>\n, -m) = ^ 5 (A^^^(iV, n, 

(C-15) 

where 

gx/iurji^: m) ■= {N + n + A, -N - m + g.\g\n + z/, -m + rf) (C-16) 

Remark C.l. The signs of the arguments of the Schur functions on the right hand side of 
mB can be reversed using (see ref. [19]) 

SA(t) = (-i)AI,;,..(-t). (C-17) 

Here |A| denotes weight of A, which in the Frobenius notation equals k+ai + - ■ ■+oik+f^i + - ■ ■+h- 

It follows from eq. ()2.15|1 that the partition function Z^r(t, n, m, t) is given by (l(l-15f) . 

where, by the Cauchy-Littlewood identity 

eEL.Er..<’4“’ = SA(t''>)s„(tPI)s„(fW)s.,(t®). (C-19) 

A,/i.,i/,77 

Using (l(J-17|l . and the product formula 

SA(t)s^(t) = ^c^^s„(t), (C-20) 
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where are the Littlewood-Richardson coefficients (see [19]), we obtain 

X,flyU,7] 

where 


h,u,{N,n,m) = (-l)l^(^+i)+-A^iV! 


''rj 7) 


'9\' 


/(iV, n, m) 


,u' ,r]' 

(C-22) 

Wick’s theorem can then be used to evaluate the right hand side of ()(hl6fl and express the 
result in terms of the bimoments. The quadruple Schur expansion follows here from the 
fact that Zjv(t, n, m, t) is essentially a r function for the two-component TL hierarchy. 

Rather than detailing this calculation, we now give an alternative method of evaluating 
I\fiur]{Xf,n,m) that uses the particular form of g in eq. ( 12 . 5 j) . and which directly yields 
the vanishing of n, m) if the length of any of partitions A, /i, z/, rj exceeds N. 

Using ()2.15j) and ()C-13|1 . we have 


= ^ SA(t(^))s^(-t(2))(iV + X,-N + /i|A(0, n, m, f)^|0, 0), 

A,/i. 

(C-23) 


where 


N 

,n,m,t)^ = n 

i=l 

It follows from (1D3) and (ID3), the sign count from interchanging the orders of fermionic 
operators, and the orthogonality relations that 





N 

{N, -N; X, n\Y[f^^\xi)P\yi)\0,0) = (-1)^^(^+^)sa([x]) v(M)^v(a:)A^(l/), (C-25) 

i=l 

which vanishes unless f'(A),£(/r*’') < N. Substituting (I(h24fl . ()(l-25fl into (l(l-23fl . using 
()(1-17|) and changing y —> in the summation then gives 


Zjv(t,n,m,t) = ^ f dg{xi,yi) 

Va=l 




w(N)s^([i/]). 


E{X),E{g.)<N 
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(C-26) 


At this point, the calculation becomes identical to the one in section |21 Using the 
Cauchy-Littlewood formula 


^V{x _ 


^ Sy{[x e 


V{y-\m) _ 


Sy{[y ^])Sr,(t^^^), (C-27) 


^{u)<N 


V 

i{r])<N 


the identities (see section El 

= {xi---XN)~'"^Si^{[x]), S^([2/"^]) = (2/1 ■ ■ ■2/iv)“''^s^([2/]) (C-28) 

where z> = (z>i, ..., z>Ar), p = (pi, ..., pat) and 

p. ;= ui - z/jv_i+i, fji := r]i - p7V-i+i (C-29) 

and the product formulae 

SA(N)si;([a;]) = Yc%Sc,{x), s^([2/])sr)([2/]) = 5^c2^jS/3([2/]), (C-30) 

a f} 


we obtain 


ZAr(t,n,m,t) = Y (C-31) 

where 

h,.,l(N, n,m)^Yl (C-32) 

a, 13 

and, in terms of the bimoments Bij, 

■= iV!det(5;,+„,;.+™)|i<,,,<^, 

li := tti — i + N — Ui, I'- := Pi — i + N — rji. (C-33) 
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